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Abstract 

Revisiting canonical integration of the classical solid near a uniform ro- 
tation, canonical action angle coordinates, hyperbolic and elliptic, are 
constructed in terms of various power series with coefficients which are 
polynomials in a variable r 2 depending on the inertia moments. Nor- 
mal forms are derived via the analysis of a relative cohomology problem 
and shown to be obtainable without the use of cllitptic integrals (un- 
like the derivation of the action-angles). Results and conjectures also 
emerge about the properties of the above polynomials and the location 
of their roots. In particular a class of polynomials with all roots on 
the unit circle arises. 
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1 Overview 



Integration of a rigid body motions is well known. In this paper we recall 
in Secj2] the kinematic description of the coordinates in which the system 
is obviously integrable by reduction to quadratures (Deprit's coordinates, 
described in Eq. (|2.5p ). We try to treat simultaneously the motions near 
all proper rotations, at the expense of having to append labels a = ±1 to 
most quantities thus making the notation a little heavier. This is done to 
keep always clear that the stable and unstable motions are in some sense 
described by the same analytic functions. 

Our first aim is to obtain canonical coordinates p, q, in the vicinity of 
the proper rotations, such that the Hamiltonian, at fixed total angular mo- 
mentum A, becomes a function of the product pq (near unstable proper 
rotations) or of p 2 + q 2 (near the stable ones). Although existence of such 
coordinates is well known and their determination is in principle controlled 
by appropriate (elliptic) integrals, [H Sec. 50], it remains, to our knowledge, 
only implicit in the literature, j2[ Sec. 69]. 

In Secj3]explicit expressions for the elliptic integrals are derived, Eq. (|3,7p . 
which give the motions in terms of variables p' , q' evolving exponentially in 
time on the hyperbolae p'q' = const or rotating uniformly on the circles 
p' 2 + q' 2 = const and "reducing" the problem to determine the nontrivial 
scaling function C (depending on p'q' or p' 2 + q' 2 ) which leads to the above 
mentioned canonical coordinates viap = p'C, q = q'C. In Appendix lAl given 
details on the computation of appropriate elliptic integrals: this should also 
clarify the advantages of the cohomological analysis of the following SecH] 
and[6j 

In SecJH the scaling function C is related to a Jacobian determinant 
between two differential forms, and it is computed: this is done, avoiding 
evaluation of other elliptic integrals, by computing the cohomology of the 
forms dB Adf3 (where B, j3 are two Deprit's coordinates, see Secj2]) and dp' A 
dq' relative to the functions p'q' or p' 2 + q' 2 . The canonical coordinates are 
therefore completely determined, Eq. ()4.10|) . (l4.1ip (together with Eq, (l3.7p ). 

The latter expressions are still quite implicit and in Secj5]we bring them 
to a form suitable for the evaluation of the integrating coordinates via com- 
putable power series; and as an example a few terms of the scaling functions 
C and of the normal form of the Hamiltonian are evaluated, Eq. (|5.4p . (I5.6p . 
Doing so various families of polynomials, in a variable r 2 depending on the 
inertia moments, arise. The polynomials appear to have zeros on very special 
locations suggesting possible conjectures. 

To understand the properties of the polynomials which arise in the con- 
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struction of the normal forms we show, Sec|6l that the normal form can be 
generated by solving another problem of relative cohomology: this leads to 
an explicit determination of the normal forms in terms of a further family 
of polynomials. In this case we prove that some of the new polynomials 
have all roots on the unit circle (with the help of theorems on the location 
of zeros of symmetric polynomials): and we conjecture a relation with the 
Lee- Yang theorem for the zeros of "ferromagnetic polynomials". 

In SecJT] we determine the action angle coordinates for the other two 
Deprit variables (which either do not appear in the Hamiltonian or are 
constants of motion) . 

In Secl8]we apply the relative cohomology method of Sec[4jand[6]to the 
pendulum, to illustrate the cohomological method. And to the geodesic flow 
on revolution ellipsoids, with the purpose of showing that in the latter case 
also appear natural families of polynomials in a parameter r 2 (depending on 
the ellipsoid equatorial and polar axes): but in this case the zeros are not 
located on the unit circle, although almost so: Eq. ([8.6p . 

Appendices contain supplements (e.g. the proof of Chen's theorem). 

2 Solid with fixed center of mass 

The theory of Jacobian elliptic functions, for reference see [3], yields a com- 
plete calculation for the motion of a solid with a fixed point. This is revisited 
here, to exhibit a few interesting properties of the relevant elliptic integrals. 

The Hamiltonian of a solid with inertia moments Ji, J2, J3, in Deprit's 
canonical coordinates (K z , A, B,j,ip, B), see problems in Sec. 4. 11 of [3], is 

1 B 2 1 ,s\n 2 B cos 2 /?, 9x 
H(K z ,A,B n ,<p,P) = -— + -(-^ + -^) (A 2 - B 2 ). (2.1) 

and it will be convenient to order the proper axes of the solid body so that 
J3 < I\ < Ii- In this way the motions with B = and 3 = 0,tt or ±^ 
are stable or unstable rotations, respectively. For convenience the unstable 
rotations will be studied by shifting by 5 the origin of the angle 8, so that 
for the unstable rotations we shall use the Hamiltonian: 

H(K„A, B, % = \ I + + **£.) (A* - B 1 ). (2.2) 

Motions can be described by referring them to the two node lines m and 
n: m being the node between the angular momentum plane, orthogonal to 
the angular momentum M = Ak, and the fixed reference plane i — j, and n 
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being the node between the angular momentum plane and the inertial plane 
1 — 2. The angle 7 locates the node m on the fixed plane and the angle <p 
locates n on the angular momentum plane with respect to m. 

k = z 



i,j,k=fixed frame 

i,j,k=fixcd angular momentum frame 
ii)i2, i3=comoving frame 




1 = x 



1 = mr n't 

The Deprit angles. Herexi is the node line (i, j) n (ii, 12), n is the node line (12,13) H 
(i,j) and m = i is the node (i,j)n(i,j). The j axis is not drawn and k || M. Colors 
are meant to suggest that the angles 7, if, f3 lie on different planes. 



With reference to the Hamiltonian in the form Eq. (|2.2p . where (3 = 0, tt 
correspond to unstable rotations around the ii proper axis, the pairs of 
canonical variables are (B,(3), (A, if) and (M z ,j); sin# = while the angle 
5 is constant. At the "fixed points" B = 0, (3 = 0, vr or B = 0, /3 = ±§ the 
angles (p and 7 rotate at constant speeds: (p = and 7 = or up = y 
and 7 = respectively; the angle between the angular momentum and the 
inertial axis 13 is ^ and 8 = 0. For the stable rotations the role of the axes ii 
and 12 are exchanged. Hence the coordinates are adapted to describe proper 
rotations around the axis 1 (/3 = 0, n, unstable) of intermediate inertia or 2 
(/3 = ±7^, stable) of largest inertia. 

The inertia moments I3 < Ii < I2 will be used to define the quantities 



J+ = h, 

t-1 def 1 1 
J 31 — J 3 

j def 
J+ — J31, 



1 5 



j def 
J- = *2, 

t-1 def 1 _ 1 
J 32 — J 3 ' 

T def T 
J- = ^32 



2 > 



J 



_1 def 1 
12 — 1 1 



(2.3) 



The Lyapunov coefficients at the unstable fixed point are ±A+ real with 
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A i = - — — r and at the stable are ±A_ imaginary with A_ = - %A — r . 

(J+Jl 2 )2 (J"-Jl 2 )2 

We try to use a notation valid for both the stable and the unstable ro- 
tations: therefore a label a = ±1, often abridged into a = ±, is appended 
to most quantities: with a = +1 referring to the unstable case and a = — 1 
referring to the stable one. Calling U the total energy, define for a = ±1: 

hi = f 2U - A2 J J \ f a m (2U - A'/Jaft 1 S b\ (J" 1 ) 2 , 
J a 

„1 d ^l Al J 12 Jjdef J 12 2 d £_f as l ~ ar l (0 A \ 

2{7 — AVJa' a ~ Jo" 1 ' a ~ 1 + asl ' 1 j 

n 2 

2 —2/1, 2\ i2/i i 2n i/a 



Or 1 ) 2 

then g a is a characteristic inverse time scale and, referring to Eq. (|2,2p for 
a = + and to Eq, (12.1j) for a = — , 



9 2 =A 2 J 12 1 J- 1 (a+ ' 



s 2 ' 



$ = ± g a V (1 + sin 2 /3)(1 + ar 2 sin 2 /3) (2.5) 
B=±b a 



1 + as 2 sin 2 /3 _ o a / 1 + as 2 sin 2 /3 



1 + ar 2 sin 2 j3 J^ 1 \ 1 + ar 2 sin 2 /3 
with .A, -ftT^ = j4 cos (5 constants of motion. If a = — it is (2U — j-) > 0. 
Suppose first (2C7 — > and small: then 

* .«2>o. «>o, {H^- 

'I » 1, { £ >*t < 1 , 1 + or2 = £| = ^, (2.6) 

2 _ fc 2 + ar 2 b\ _ ggJg _ 2 ~ ~ 

a " 1-A: 2 ' l-fc 2 "l + ar 2 -^ J+J - 

and motions near 1? = 0, /3 = and with energy U close to A 2 /2J a and can 
be expressed in terms of Jacobian elliptic integrals (see [3] for notations). 

Remark: Rotating the positions of the inertia moments in H the meaning 
of the variables changes but the following analysis remains essentially un- 
changed (up to renaming variables): it appears in this way that the proper 
rotations around axis ii are (linearly) unstable while the rotations around 
the highest and lowest inertia axes (i 2 and 13) are (linearly) stable. 
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3 Hyperbolic and elliptic coordinates 



"Linearization" of rigid body motions is, of course, well known to lead to 
elliptic integrals, see [21 p. 144]. Here we rederive the integration in a form 
suitable for our purposes of action-angle coordinates construction. 

Let \J— 1 = — i so that g a = y/a\g a \; for a = + Eq. (|2.3p describe the sep- 
aratrix branches emerging from the proper rotation around ii controlled by 
Eq. (|2.2p . while for a = — they describe the motion near the stable rotation 
around i 2 controlled by Eq. (|2.ip , 

Let U — A 2 /2J a > close to 0; from the definitions of the Jacobian 
elliptic integrals, see [3 2.616], and setting 



K = \A ~ kl = y I * , u a = sjl + as 2 a g a t = g a t (3.1) 

it follows (choosing the sign + in the Eq. (|2,5p . for instance) from ([3j 8.153.3] 
for a = + and [3,, 8.153.1-8] for a = — after correcting the obvious typo in 

[3j 8.153.8] and by k'_ = —i^J k 2 _ — 1), see Appendix Rl for details, 

Bit) = ba = b cn H^X) 

dn(u a ,k a ) a dn(-iu a ,k' a ) , 

2 _sn 2 (-m a ,fc^) 

gin2 m = sn { u k a ) = B| ggg (3 . 2) 

a a cn^(— lU a ,k' a ) 

■ _ (1 + asl)g a dn(u a ,k a ) __ dn(-iu a , k' a )cn(-iu a , k' a ) 



l + aslcn 2 (u a ,k a ) a 1 - T ±^ S n 2 (-iu a , k' a ) 

where the last relation is deduced from the equations of motion and the 
expression for sin 2 j3. 

Notice that 0(0) = and if (2U - A 2 /J a ) -> 0+ it is -»■ 1 T , -»■ 
J^ 1 ^ 1 ^ 2 , 6 a - >■ + and motions with this energy are "like" the motions 
close to the separatrix or, respectively, close to equilibrium of a pendulum. 

For (2U - A 2 /J a ) < 0, or for k a > 1 hence k' a imaginary, the above 
relations have to be interpreted via suitable analytic continuations. Some 
of the following formulae become singular as U — > A 2 /2J a : the singularity 
is only apparent and it will disappear from all relevant formulae derived or 
used in the following. 

Introduce variables useful in the following (most of them appear in the 
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theory of Jacobi's elliptic functions): 

7T 

K(k) d = f [ 2 (l-k 2 sm 2 9)~he, 
Jo 



x' a d = f ae-^, I 2d = f 4J + L ( 3 - 3 ) 

2 ;/2 — 2/i , 2\ j. I 1 7F 



2K(k' a ) 



The g a ,go, a depend on k a and can, and will, be imagined as functions of x' a . 

Let a a d = arcsin , 1 = —iv a = [n a , de and the integral can be 
evaluated by series, [31 2.511.2] and [31 8.113.1], defining P m and V 

- iVa =?^a a -^if: H)(-k' 2 rPm(—^- 2 ) 

m=l \ / 1 a (3 4) 

4f«(«.+ /^( 7 f i ,tf)) 

7T V v 1 + ar~ / 

with P m (z) = ^(EfcU 1 2^^=^p7! This is wel1 defined near tne 

proper rotations , where k a — > 1, as /c^ tends to 0. Therefore, as derived in 
detail in Appendix lAl 

la = T(arlX 2 )= f e^ = (i + ^) e ^(^)^^ 

Vl + ar^ Va 

and |7 + | = 1, 7_ = ie x with A real and j_ real, with A real. 
Introduce the new variables p' + ,p'_ defined by 

„' [Zl p a\/ago, a t 



p' a =d y Jx' a e a " /a90 < at , a = ±, hence x' a = p' + p'_ . (3.6) 

The above relations can be algebraically elaborated into the transformation 
of coordinates (B,(3)< — >(p' + ,p'_), see again details in Appendix [Al 

d f 00 a m x' m p' 

B=R( P ' + , P '_) = alg , a (x' a ) E E 1 + a( ;™w )2 

fj,=±l m=0 v ° ^f*' 

P=S(j/ + M= f ajT l E ^« n arctanh(x' a >^) 

n=0 (T,7j=±l 

and in the new coordinates the flow t —> p' a e'^°' a90 - at generates a solution of 
the equations of motion if \p±\ is small and if the real parts of p'± are > 0. 
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However the equations of motion are analytic and so are the Eq. (j3.7j) : hence 
Eq. (|3.7p yield a solution under the only condition that \p±\ are small, to be 
required for its convergence. 

Eq. (|3.7p give a complete parametrization of the motions but the coordi- 
nates Pu are not canonical: the determination of the canonical coordinates 
will be discussed and done in SecUl 

Remarks: (1) via identities relating Jacobian elliptic functions it would be 
possible to avoid considering functions that are defined by analytic continu- 
ations, for instance when the modulus k 2 . > 1 or k' 2 < 0, and reduce instead 
to considering only cases with the "simple" arguments (expressing functions 
by elliptic functions of other arguments, like h 2 = 1 — k~ 2 , h' 2 = k~ 2 , when 
necessary): we avoid this because it would hide the nice property that the 
stable and unstable cases can be seen as analytic continuations of each other. 

(2) The formula for (3 reminds of one found by Jacobi which he commented 
by saying that " inter formulas elegantissimas censeri debet", [5, p. 509] (i.e. 
" it should be counted among the most elegant formulae ", see also [6]. 

4 Canonical coordinates construction. Jacobian 

The motions p' a — V p' a e CT ^ 9o - at solve the equations of motion if the real and 
imaginary parts of p'± are positive. But the equations of motion are analytic, 
hence the formulae of the previous sections give solutions of the equations 
independently of the sign of p' a , provided the series converge. Convergence 
requires \p'-\ <C 1: which represents many data, in particular those in 
the vicinity of the separatrix or of the stable proper rotation. 

The coordinates can be called "hyperbolic" or "elliptic". We also see 
that time evolution preserves both volume elements dBd/3 and dp'dq'; which 
means that the Jacobian determinant f,^'^. must be a function constant 

over the trajectories, hence a function D a (x') of x' d = ' p' + p'_ in the two cases. 
Note that D a (x') has dimension of an action. 

Introduce the variables p' , q', related to p'± in Eq. (|3.6p by: 



to deal with real coordinates when useful; it is then possible to change coordi- 
nates setting p± = C a (x' a )p±, (or p d = C a (x' a ) p' , q = C a (x' a ) q') and choose 
C a (x') so that the Jacobian determinant D a (x') for (B,(3)< — > (p+,p~) is 
= 1. A brief calculation shows that this is achieved by fixing 




(4.1) 



S 



CaV) = \, f D a (y)dy, (4.2) 

X J o 

which is possible for x' small because, from the equations of motion it is 
D a {0) finite (e.g D + (0) = 4/5(0) sin > 0). Therefore map (B,/3)<—> 

(p,q) is area preserving, hence canonical The Hamiltonian Eq. (|2.ip be- 
comes a function U a {x a ) of x a = and the derivative of the energy with 
respect to x a has to be go, a (x' a ) (because the p+,P- are canonically conju- 
gated to B, 0). Note that x a has the dimension of an action, while p, q are, 
dimensionally, square roots of action. 

This allows us to find D a (x'): by imposing that the equations of motion 
for the (p, q) canonical variables have to be the Hamilton's equations with 

Hamiltonian U a {x a ) d = U a {x' a ) = H{B,p) it follows that ^ga) = g , a {x' a ), 

A p dU a (x' a ) dx' a _ dU a (x' a ) _ _d_( T t n ( T l ^No- 1 ^ 2 _ r> ( T > \ 

l - e - dx' a - 90,a or — - go,a d^r{X a O a {X a ) )Z 2 _ g Qa U a (x a ) 

by the above expression for C a {x' a ). The just obtained relation gives 

D a(x' a ) = goA^a)' 1 -^rU a (x' a ) (4.3) 

which is an expression for the Jacobian g^fy — d(^',q^) = Wrftf) ( no ^ e that 
the Jacobian between (B,/3) and (p,q) is identically 1 by construction). 

To proceed it is necessary to determine the function C a (x' a ) 2 . The idea 
in [6] is to make use of the general theory of relative cohomology classes to 
determine D a {x' a ) and, therefore, to find a complete expression of the action 
angle variables following a procedure employed in the theory of limit cycles 
of planar vector fields, 0[5]. 

A derivation can be based directly on the equalities on symplectic forms: 

dBAd!3 = D a (x' a )dp' A dq' = J-[ x ' a C a {x' a f] dp' A dq' 

dx' 

A (44) 



() 



2 ; 2 -j^\x a Ca{x a f}dp + f\dp'_ 



'-a 



To simplify the notation we drop for a while the label a from x h s. The idea 
is to compute only the cohomology class of the volume forms in the relative 
cohomology of the function x' = p' + p'_ , |8] . 

For any 2-form (f(p' , q')dp' Adq' with ip(p' , q') = J2 m n fm,nP /m q' n the "co- 
de/ 

homology class" relative to x' = p' + p'_ is the function defined by ■ip(x') = 
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J2m fm,m(p'q') m - Therefore an expression for the Jacobian D a (x') is ob- 
tained by finding an expression for the cohomology class of the form dB A d/3 
relative to x' = p' + p'_ . 

Remark that Eq. (|3.7p imply 

+00 i±2 

R'(p'+,p-)dS'(p+,p'-) = igo(x') E E — 1 — 

m,n=0 n,o,ri=±l (4.5) 
m+n /m /n-r, 

x 7 « v dv' +0 

(1 + a(x"X)2) (1 _ (a/XTS) 2 ) M a 

where is a form with cohomology (i.e the first term is RdS "modulo 
d(p' + p'_y j because the cohomology of d\G(p' + p'_)d(p' + f/_) \ is, for all G's, 0). 

The terms with p + dp + and p_dp_ do not contribute to the cohomology 
and, up to terms (denoted 0) not contributing to the cohomology, the form 
in Eq. (|4.4p becomes (as only terms with [i = —a contribute) 

R'(p' + ,p'_)dS'(p' + .p'_) = I 9 o(x')- £ E *' G<hL 



* m,n=0 r],a=±l P* 7 (4 6) 



(1 + a(x' m p'_ a ) 2 ) (1 - (x /n 7VJ 2 ) 

If RdS = x>F(p> + ,p'_)\( d ^- - |±) + and F(p' + ,p'_) = Eh,k=o Fh,kP'lp*, 
only the terms with k = h contribute and we can take 

1 00 a m+n x' m+n+l ri^^f 

F = **^i^JL (1 + a(*"X) 2 ) (1 - {x' n llp^Y) (47) 

hence the cohomology with respect to p+p- is that of (^- — times 
2I 9o( x ') V" x'r]^a m+n x' m+n j v (x' m+n+1 7 v ) 2h ^^ 

^ r?=±l m,n,h 

rfr_ ~ (/ + i) 

which contributes 



(4.8) 



i go (x') *'E— E r+a(x , ( ^ )2 
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because the cohomology oi d(x' f (x')(^- — ^-)) is —2-^j(x'f(x'))dp' + Adp'_. 



l-a 



Therefore, by Eq.flOJ and by dp' + A dp'_ = (f)~dp' A dq', see Eq. (j3~3j) : 
C a (*) =2/ 50 , a (x) £ (-) L^-pTi)^ 

r;=±l £=0 1 + ^ 7«J 

= 9f f*\yrRyr (^ + iy^r(org,^)" (4.10) 

Summarizing we set, with defined in Eq. (|4.ip . (|3.6p . 



r a)> 

/ / . P+P- i rat i 2n ( 4n ) 

s =/(p+p_), 2^T =a: « = aJ a C (ax a ,ar ) 



where p + p'_ = pV if a = + and p' + p'_ = (p + q ) if a 



5 Evaluation of Jacobians and normal forms: 

It it possible to evaluate C 2 and the normal form for the energy U = H 
in the canonical coordinates (p+,P-) up to large orders in x a . This is done 
simply by deriving power series expansions, in x' a first and then x' a in terms 
of x a , for the functions in Eq. (15.31) . 

Consistently with Ea. ([4"TTj) C 2 can be expanded into a power series of x 
or x 1 whose coefficients are functions of the parameter r 2 . There is a simple 
relation between the cases a = ± and the coefficients of the expansions in x 
can be computed exactly with a finite computational algorithm up to any 
prefixed order: this is described below. 

The analysis of the Jacobian C a (x,r) 2 can be based on the computable 
expressions of 7 y a in Eq. (|3.5p . 



Normal forms express H a = 2H ^_( a . From Eq. (|2.ip , (|2.2p we derive 

aJ 12 

A 2 s~ 2 (see Eq. (|2.6|) ) in terms of k 2 and, by Eq. (|3.3p . x' a in terms of k 2 : 

.K(fcq) 



ar„ =ae K ^) , a = ± i 5 - 1 ) 



2 V 

which can be inverted in terms of the functions F(z) = -^=^ 

F(ax')\^def 



J2„-1) 



2 



tt as 



1+2 V °° nZ (2n)^ 

^ ( ^ )2 = (lT^) 4 = V(ax,) ' a = ±1 (5 ' 2) 



li 



as shown in [3, 8.198]. Therefore the above relations give explicitly the 
functions H a in terms of x' , hence they determine the functions U a (x'). 

To express H a in terms of x we have to determine the Jacobian function 
and invert the relation x = x'C^(x') in the two cases. This is done by 



where the relations r 2 _ a = pf^p ; ~jir = ^ r °' following from Eq. (|2.6p or 
(|2.3p . have been used and the function V is defined in Eq. (|3.5|) . 

This is remarkable because the summation over rj = ± is analytic in ar 2 
in spite of the fact that each addend has a branch point in y 'ar 2 at r a = 0. 
Therefore Eq. (|5.3j) shows that the two cases a = ± are described by the 
same functions (of course evaluated at completely different points) and the 
power series for C and U can be computed setting a = + and then the case 
a = — is obtained by replacing x by — x and r 2 by — r 2 , ^4r + by Ar_. The 
simple relation between the unstable and stable cases implies that that we 
only need to study one of the two cases: we shall only make computations 
in the hyperbolic case a = +. 

At this point the Jacobian and the normal forms can be computed as 
formal series in x after solving the implicit function problem x = x'C 2 (x' , r 2 ) 
in the form x' = xC (x, r 2 ) (in a power series in x). 

An algebraic numerical evaluation for a = + (unstable case with, for 
simplicity, r = r + ,x' + = x') yields results that for n < 7, if £ = ^ , are 



remarking that go,a = V° #a 2K^fc / ) ana - U a {x') can be written as 




(5.3) 



.2 



a; C(x, r 2 ) 2 = 4£ - 24(-l + r 2 )£ 2 + 32(3 - 14r 2 + 3r 4 )f 
-16(-1 + r 2 )(19 - 242r 2 + 19r 4 )£ 4 
+24(35 - 1140r 2 + 3026r 4 - 1140r 6 + 35r 8 )£ 5 
-192(-1 + r 2 )(ll - 740r 2 + 3426r 4 - 740r 6 + llr 8 )£ 6 
+64(77 - 10206r 2 + 103635r 4 - 211460r 6 + 103635r 8 
-10206r 10 + 77r 12 )£ 7 + ... 



(5.4) 



12 



Finally the polynomials in r 2 in Eq. (j5.4|) . coefficients of order n > 1 for the 
power expansion in £ n , can be conjectured to have all roots real: this has 
been numerically checked for n < 15. 

The normal forms of the Hamiltonians (i.e. their expression in the vari- 
ables p,q) are derived by expressing s 2 in Eq. (|2.ip in terms of r 2 and k 2 

and then expanding k in powers of x' using x = e K (* ) in Eq. (|5.ip , via [H 
8.198] (notice that in our hyperbolic case x' is the q of the quoted tables with 
k and k' exchanged and finally computing x' = xC (x,r 2 ) by inverting the 
Jacobian relation x = x' C 2 (x',r 2 ). This can be implemented algebraically 
on a computer. 

With the definitions in Eq.flO]) H a d = 2H ~f_{^ , in the unstable case 

a J 12 

(a = +1) or in the stable (a = —1) case, can be expressed as function of 
x = 2 Ar+ or ' respectively, of x = p 2 ^ r q obtained by solving x = x'C 2 (x',r 2 ) 
as x 1 = xH(x,r 2 ) and replacing x' by S (x , r 2 ) in U (x 1 ) in Eq. (|5.3p (notice that 
J^ 1 = r 2 Ja)- The resulting functions are normal forms for the Hamiltonian 
since p, q are conjugated variables. Hence 

H(x a ) = ^T + n ( aX ^ ar a) ( 5 - 5 ) 

for a = ±1 where H(x,r 2 ) d = f ^ = , where W(x) is in 



Eq. fl5.2|> . One obtains 

H(x, r 2 ) =4x- 2(-l + r 2 )x 2 - (1 + r 2 ) 2 x 3 
_!(_! + r 2)(! + r 2)2 x 4 - 1(1 + r 2 ) 2 (ll - 10r 2 + llr 4 )x 5 
7 

- — (-1 + r 2 )(l + r 2 ) 2 (3 - 4r + 3r 2 )(3 + 4r + 3r 2 )x 6 (5.6) 

— -(1 + r 2 ) 2 (527 - 332r 2 + 330r 4 - 332r 6 + 527r 8 )x 7 
64 

9 

(-1 + r 2 )(l + r 2 ) 2 (1043 + 548r 2 + 1058r 4 

512 

+548r 6 + 1043r 8 )a; 8 + . . . 

It is remarkable that the polynomials in r appearing as coefficients of order 
2 < n < 15 in the power series expansion in x of the normal form have, for 
all n 's, just roots of unity as it appears by evaluating them numerically. 
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6 Normal form without determination of action- 
angle variables 

The function T~L(x) being the normal form of the Hamiltonian in action- 
angle variables is unique, while the action-angle variables themselves are 
not unique. Therefore it makes sense to see if the normal form itself could 
be derived without actually determining the angles (notice that the actions 
being a function of the energy are uniquely determined). Here we show 
that this can be achieved quite directly without having to deal with elliptic 
integrals: while in the elliptic case, for instance, the classical prescription, 
[TJ Sec. 50], would be to study the Hamitonian as a function of the integral 
§B& with B defined in Eq.flZS])). 

As it is well known the full determination of the action-angle coordinates 
is much more difficult than the determination of the normal form of the 
Hamiltonian, see for instance the "conjecture" in [9 J and its proof in [10] . 
Therefore a direct determination of the normal form is worth investigating: 
and as it will appear here, the construction yields unexpected insights into 
our problem. 

Consider first the normal form near the stable rotation. In contrast with 
our first method, we do not start with the Jacobi's coordinates in which the 
motion is linear. 

The Hamiltonian value U will be written as (see Sec]21 a = — ): 

2U-A*/h = B 2 {1 _ r 2 gin 2 p) + ^2 gin 2 p 

For simplicity of notation the label a = — will be omitted. Introducing X = 
B^J\ — r 2 (sin /3) 2 , Y = r A sin j3 the Hamiltonian becomes U = ^j^i^ 2 + 

Y 2 ) + ^ and the form oj = dB A d(3 is, if X± d = ^f, 

1 dX A dY 



dB Ad/3 



(l-r- 2 (^) 2 )(l"(^) 2 ) Ar 




r> / 2n (X + X-\ndX A dY n , _ . 2 , 



n=0 
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and is a form with 0-cohomology relative to X + X^. 

Therefore, see proposition l].|12j. or Theorems 1.1, 1.2 in |13| . there 
is a change of coordinates X, Y — > X , Y which, if H = (2Ar)' i ano - 
H = (2Ar)" 2 i ^ ee P s H e< l ua l to H' and makes uj proportional to dX' A dY'\ 



oo 



n=0 



n / 2n rr,n dX' A dY' def _., TrA dX' A dy' 
u, = ^ Pn(r 2 )ff m ^ = J £>(H0 ^ (6.3) 



Let p = -S=d(H'), q = -¥L=C{H') with C chosen, see also SecH so that 

x 2 , „dC(z) 2 



00 p / 2\ 

^(^) 2 = E rrr^ n < ( 6 - 4 ) 



dB A d/3 = dp A dq, i.e. as C{z) 2 + ^^f^ = or 

" 2 ) 

n=0 (^+ T )' 

For instance, up to 0(z 7 ): 

z C{zf = z-\{-\- r 2 )z 2 + 1(3 + 2r 2 + 3r 4 )z 3 

_ A(_i_ r 2)( 5 _ 2r 2 + 5 r 4 ) z 4 
7 

+ —(35 + 20r 2 + 18r 4 + 20r 6 + 35r 8 )z 5 ( 6 5) 

21 

(-1 - r 2 )(63 - 28r 2 + 58r 4 - 28r 6 + 63r 8 )z 6 

128 

33 

+ (231 + 126r 2 + 105r 4 + 100r 6 + 105r 8 + 126r 10 + 231r 12 )z 7 

256 

and a normal form of the Hamiltonian is 

A 2 A 2 r 2 

U{M)= 2h + 2J^ n{xy) (6 - 6) 
where 7i(x,r 2 ) is 4 times the inverse function to x = zC(z) 2 and x = 

— 2 — 2 

p 2 ~aI ( w hh r = r_): this gives a normal form, because the variables p, q 
are canonically conjugated to B,f3. For instance (and in agreement with 
Eq.([5SD): 

H(x) =4x- 2(1 + r 2 )x 2 - (-1 + r 2 ) 2 x 3 - 5(i+r 2 )(-i+r 2 ) 2 x 4 

_ 3 (-l+r 2 ) 2 ni - n 2 ,11,4^5 _ 7(-l+r 2 ) 2 (l+r 2 ) , Q _ „ 2 , Q^4 W 6 



L(j -(11 + 10r 2 + llr 4 )x 5 - 'l-^; u+r ^ Q _ 2? , 2 + 
— -(527 + 332r 2 + 330r 4 + 332r 6 + 527r 8 )x 



(6.7) 

( ) ftm j- 332r 2 + 330r 4 + 332r 6 + 527r 8 )x 7 

- , ., (1043 - 548r 2 + 1058r 4 - 548r 6 + 1043r 8 )x 8 + . . . 



9(l+r 2 )(-l+r 2 ) 2 , 
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It should be noted that the above analysis does not determine the canonical 
map into action-angle variables of the initial coordinates and, in particular, 
a further change of coordinates which transforms the plane p, q into itself by 
rigidly rotating all circles centered at the origin by a radius dependent (and 
even r dependent) angle will lead to a normal form. For instance calling 
Q = p + iq the change of variables CJ = Qe 1 ^®^ ) is a canonical map, for any 
choice of Q(r 2 ), and it keeps the system in normal form. 

Since the Hamiltonian is uniquely cast in normal form, the polynomials 
P n (r 2 ), coefficients of the expansion of % in powers of x, have an intrinsic 
significance. It is therefore interesting that they seem to enjoy the property 
of having all zeros on the unit circle: we have checked this for 1 < n < 15: 
are really the zeros on the unit circle for all n ?. 

The location on the unit circle of the roots of our polynomials is certainly 
related to the property that the polynomials 

appearing in the above expression for C 2 have all the roots on the unit circle. 

The latter is a special case of the general result on polynomials, that 
we call "polynomials with symmetric, positive and monotonic coefficients" , 
of the form Q(z) = Eh=o F hZ h , with F h = F n _ h and F h „ x > F h > for 
< h < [2±2]: 

Theorem: (Chen, 1995) Polynomials with symmetric, positive and mono- 
tonic coefficients have all roots on the unit circle. 

The theorem can be found in [14) . and it is described, for completeness 
in Appendix [Bj 

This is analogous to the conjecture formulated after Eq. ([5.4j) . and it is an 
algebraic property which could, we think, also be seen as follows. Consider 
the polynomials 

Q^) = z n e^ Jiii ^ )2 z^i ai (6.9) 

<Tl,...,cr n =±l 

with Jij > 0. Then the coefficients Qh of z 2h in Q have the form of a sum of 
(?) exponentials of sums of some of the Jij : the sums define quantities ah = 

a n -h d = e p ( h \ By adjusting the constants Jij the parameters p(h — 1) —p(h) 
for < h < I 1 ^] can independently take any value > so that e p ^~ p ^ = 
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can be achieved by suitable choices of the Jij > 0. Therefore Q would 
have all zeros in the unit circle by Lee- Yang's theorem, |15j , achieving a 
alternative proof of the above theorem. 

The latter is a conjecture that can be checked rigorously for n < 5 and we 
found some preliminary numerical evidence that it should also work at least 
for n < 8. The proof for n = 4 emerged from a discussion with A. Giuliani; 
he also has a proof solving the case n = 5 and gave us the reference to 
Chen's theorem. 

The polynomials in r 2 appearing as coefficients of x k in the normal form 
Eq. (j6.7p are sums of products of polynomials P nj of degrees which add up 
to 2(k — 1) = Ylj n f an expression for the inverse function to zC(z) 2 can 
be found in terms of a "tree expansion" (quickly sketched in Appendix IC|) . 

Sums of products of symmetric polynomials with symmetrically decreas- 
ing positive coefficients is a polynomial with the same symmetry but it will 
not be, in general, decreasing. Hence the above theorem on the zeros of 
polynomials does not allow us any conclusion on the location of the zeros 
of the coefficients of the normal form: however we see empirically that the 
zeros are apparently still on the unit circle (up to n = 15 at least). 

This is a special property of our polynomials P n because, in general, the 
inverse function to the map x = z + J2^=i P n {w)z n+1 with P n (w) a degree 
n symmetric decreasing polynomial with positive coefficients (hence with all 
zeros on the unit circle) can be written z = x + J2^=i Qn(w)x n+1 but the 
Q n do not have, in general, all zeros on the unit circle. 

The analysis of the hyperbolic case is essentially identical: introducing 



the variable X± = B\ \ + r 2 sin (3 ± Ar sin f3 the form dB A d/3 is 



£ 7047^ (6-10) 



2W(l-sin 2 /3)(l + r 2 sin 2 /?) 2Ar t'o V ( 2Ar 



and the normal form is related to the elliptic case and given by U(p, q) = 
C ~ IDS? n (~ x ' ~ r2 )- where x = 2M and r = r +- 

Remarks: (1) The normal form in this section coincides with the one in 
SecJS] because of the normal form uniqueness. The agreement continues up 
to order 15: however the above analysis shows that this must hold to all 
orders: i.e. tths coincides with x a and our construction in SecEE] deter- 

a 

mines explicitly one among the changes (which are not unique) of variables 
(X, Y)< — >(X',Y') realizing the cohomological equivalence of the symplectic 
forms dX' A dY' and dX A dY . 
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(2) It would be interesting to find a mechanical interpretation of the prop- 
erties of the roots of the polynomials P n . 

7 Action-angles for the "momentum" degrees of 
freedom (A, cp) 

Consider the elliptic motions, i.e. rotations close to the proper rotation 
about the axis 12- Here we compute the second pair of canonical variables 
corresponding to the second degree of freedom and described canonically 

by (A, ip). Setting sn' d = sn (— iuh'^ 1 , h), cn' = cn (— iuh'^ 1 , h), dn' = 
dn (—iuh , h), 
Eq, (|3.2p implies 

. 2 a = , /2 sn /2 /dn /2 _ _1_ sd/ 2 

Sm P l-s 2 dn' 2 /cn' 2 A; 2 (l-/i 2 sn' 2 )(l- S 2 dn' 2 /cn' 2 ) A ' ' 

1 sn' 2 M-sn' 2 ) 



k 2 (1 - h 2 sn ,2 )(l - sn' 2 - s 2 (l - s 2 sn' 2 )) 

1 sn' 2 (l-sn' 2 ) _ 1 sn /2 1 sn' 2 

~ F(l-/i 2 sn' 2 )(l-s 2 )(l- sn' 2 ) ~ k 2 {l -s 2 ) dn' 2 ~ r 2 - s 2 dn /2 

and this yields 

v = r 2 + r/i " 1)sin2/3 (7 - 2) 

therefore, [3l 8.146.11], the expression of cp(t) and of the uniformly rotating 
angle ijj(t) with angular velocity ujq(A,x) (with x = x(x') 2 x') are 

^_ ( ^ )2 _2_ £ £ ( _ ir+m+ i 

rJn+m—\ 

. fjfjl i(a(2n-l)+a' (2m-l))go, a t In 

aU ( 1+x /2n-l)( 1+x /2m-l) e U 

^■ I) = ^l 1+ 1-r 2 ^ g (1 + ^-1)2 
with n,m > 1. 

Since the generating function of the canonical map integrating the mo- 
tions has the form S(A',x,/3) = A'(p + S (A' ,x, (3), |4 S Problem 4.11.5], it is 
A' = A and the angle conjugated to A must differ by a constant from ip and 
therefore it can be taken equal to ■0, i.e: 
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o 1 — r — — 

(7,ff'=±ln(T+mffY0 
v 2n-1^2m-l 



(7.4) 



OCT 



(1 + x /2 "- 1 )(l + x' 2m - 1 ) 

with z,j = (p' + icrg'), x' = p' 2 + g' 2 and = pa(x') , q' = ga(x') . 

This completes the construction of the elliptic canonical coordinates con- 
jugated to (B, /3,A,(p). The hyperbolic case, i.e. rotations close to the 
proper rotation around ii, is treated in the same way. 

The case of the stable rotations close to the stable proper rotation around 
i3 is also interesting and could be treated in a similar way leading to a com- 
plete parametrization of the motion along the unstable manifolds emerging 
from the rotations around the axis i3 which correspond A = B. How- 
ever the coordinates (B, (3) have a singularity on such rotations (notice that 
the j3 coordinates disappear from H): hence new (B,j3) coordinates have 
to be introduced. The simplest is to use the symmetry on the permuta- 
tions of the inertia moments and make use of the equivalent Hamiltonian 
H(B,f3) = \% + (^ + ^)(A 2 - B 2 ) in which the coordinates B,P 
have a suitably different meaning but can be used to describe the motions 
close to proper rotations around the axes i3 and ii (respectively stable and 
unstable). The analysis is then a repetition of the case a = + with the 
Hamiltonian Eq.(E2). 



8 Pendulum and ellipsoid of revolution 

The method of Secj6] can be applied to the case of the pendulum {H = 
+ Ig 2 (l — cos/3), [6]) or to the geodesic motion on a revolution ellipsoid 

( H = 2(fe 2 sin" 6+a 2 cos 2 9) + 2a4in 2 B ' K (4.12.6)]). 

The first case is considered in an attempt at proving the sign property 
of the normal form coefficients suggested in [BJ. In the second case it is of 
interest because the normal form has the form of a series with coefficients 
polynomials in the ratio r 2 = ^ pf the equatorial axis a to the polar axis 
b and in this case the zeros of the polynomials do not seem to have special 
properties (in particular they are not located on the unit circle). 

The normal form at the stable equilibrium of the pendulum is easily 
evaluated (details are skipped) to be 8Ig 2 'H( p 1 ^ I q g ) if %{x) is the inverse 
function to x = zC(z) 2 with: 
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oo 

C\z) = f 




4 n (n+ 1) 



(which is F{\, ^;2;4z), see [3 9.100]). And to order 0(x 8 ) the U is: 

„„ N 1 2 1 o 5 4 33 . 63 fi 527 7 9387 8 

■H(x) = x - -x 2 - -x 3 x 4 x 5 x 6 x 7 x 8 - ... (8.2) 

y ' 2 4 16 64 64 256 2048 v ; 

This is consistent with the first few terms computed in [6, (C.8)]: due to dif- 
ferent meaning of the symbols the relation between the function %{x) above 
and W(x) in the quoted reference is W(x) = ^H{Ax). However constructing 
%{x) via the algorithm in Appendix [C] the sign property suggested in [6] 
that Tiix) might have all coefficients, but the first, negative does not follow; 
see the final comments in Appendix O 

In the case of the geodesic motion near the equatorial circle of the ellip- 
soid let a, b be the equatorial and polar radii of an ellipsoid of revolution. 
The Hamiltonian for the geodesic motion is 

1 B 2 1 A 2 

~ 2b 2 cos 2 P + a 2 sm 2 /3 + 2 a 2 cos 2 (3 (§ ' 3) 
where A is the polar component of the angular momentum and (3 = 9 — 
| the polar angle, [H Problem 4.11.5]. Let r 2d = £ d = r 2 - 1, X = 
, B „ ,y = -^tan/3, X± = X±iY; then H = \{X 2 + Y 2 ) + ^ and: 



ab \ 1 + esm z ft 
A (l + tan 2 /3)§ A (i + (^)2)2 



3 



A 

a 




(8.4) 



A; 



2 \r 2h dXAdY + 



Therefore the normal form is H = ^ + 2 -^-T-L{x) with x = 4^5 (p 2 + <7 2 ) and 
%(x) the inverse function to x = zC 2 (z) with C 2 (z) = J2™ =0 P n {r 2 )z n and 
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which is interesting as it shows that the location of the zeros on the unit 
circle is a peculiarity of the particular case of the rigid body. 

Nevertheless a numerical investigation for n < 1000 with an algebraic 
manipulator indicates, as a reasonable conjecture, that the zeros become 
asymptotically located on the unit circle: for n large they are within an 
annulus of radii r' 2 < r" 2 with: 

\' J l-r°°\ + A • ^1,^.88 (8.6) 
with log a" ~ 1.992, c" ~ 0.880, log a' ~ 0.462, d ~ 1.000. 



A Hyperbolic and elliptic coordinates: details 

In this appendix we give details for the derivation of the formulae in Sec. El 
Let = —i so that g a = ^Ja\g a \; for a = + Eq. (|2.3p describe the 

separatrix branches emerging from the proper rotation around the axis ii 
and for a = — they describe the motion near the stable rotation around 
the axis i 2 . For U — A 2 /2J a > close to and from the definitions of the 
Jacobi elliptic integrals, see [5J 2.616], it follows (choosing the sign + in the 
Eq. (|2.5p . for instance): 

■ - - Z nir , Z / 



'l + as 2 a sin 2 f3(t) + as 2 a cn 2 (g a t, k a ) 

(A.l) 

and define h a , h' a in terms of k a , k' a so that k + ,k' + , /?._, h'_ G (0, 1) as: 

Then (O 8.153.3] for a = + and [3J 8.153.1-8] for a = — after correcting the 
obvious typo in [3, 8.153.8] and by k'_ = —i^J k 2 _ — 1): 

d/ja b a , cn(— iu a , k a ) . . u i—\ , \ 

B{t) = dn(u k ) = ba dn( -in k' ) = h ^~ m - h a > h *)> 

2 _sn 2 (-m a ,fc^) 

gin2 m = sn (u a k a ) = c ^gg 

1 + asf;cn z lua, k a ) 1 + as;j - 2 , — ttt 

a v **' o, cn^(— iu a ,k' a ) 
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■ _ (1 + as 2 a )g a dn(u a , k a ) __ dn(-iu a , k' a )cn(-iu a , k' a ) 

P\ T ) — 1 , „„9„_9/„. l \ —9a' 



l + aslcn 2 (u a ,k a ) a 1 — jj^?sn 2 (—iu a , k' a ) 

where the last relation is derived from the equations of motion and the pre- 
vious expression for sin 2 /3. Also the second and third relations in Eq. (|A.3P 
could be given expressions involving the moduli h a ,h' a so that one could 
choose to use the formulae with k a for a = + and with h a for a = — . 

Remark also that, [2 8.128.1], 8.128.3], K(k') = h'K(h), and also that 

h'Jk'J _ 1 
h- ~ L - 

Notice that /3(0) = and if (2U - A 2 /J a ) -»• 0+ it is -> 1^, ag 2 a -> 
J^ 1 J^2 A 2 , b a — > + and motions with this energy are "like" the motions 
close to the separatrix or, respectively, close to equilibrium of a pendulum. 

For (2U — A 2 /J a ) < 0, or for k a > 1 hence k' a imaginary, all the above 
relations have to be interpreted via suitable analytic continuations. Some of 
the following formulae become singular as U — > A 2 /2J a , but the singularity 
is only apparent and it will disappear from all relevant formulae derived or 
used in the following. 

Denoting cn = cn(u a ,k a ), cn' = cn(— iu a , k' a ), dn' = dn(— iu a , k' a ), sn /d = 
sn(—iu a , k' a ), from the equation of motion (and dn 2 - = 1 — A; 2 sn 2 -) it follows 

9a sn(-iv,k' a )dn(-iu a ,k' a )cn(-iu,k' a ) 
sn(-iv,k' a ) 1- j^sn 2 (-iu a ,k' a ) ^ 

identically for all v: convenient choice of v is v = v a with k' 2 sn 2 (— iv , k' a ) = 
(1 + OS 2 )" 1 , i.e. m(-iv a ,k' a ) = and -iv a d = am(arcsin ,K)> 

then (notice that g 2 k' 2 =g a (l + ar 2 )): 

P= ^5a(l + ar 2 )5 V sn(-»(w +r)V a ),k' a ) (a.5) 

Z TJ=±1 

which follows by the addition formulae, [3j 8.156.1], from the last of Eq. (|2.5p . 

If a a ^arcsin . 1 . i.e. e fa " = it is -»«„ = f" a , de , 

and the integral can be evaluated by series, [5J 2.511.2] and [3J 8.113.1], 
leading to Eq. (|3.4[) . (j3.5|) . As remarked in Sec|3] 1 = 1, 7_ = ie x with A 
real and 7_ real, with A real. 



It is useful to check that if I = y 4J + J_: 
50 ' a ZK(A4)' «A4) 'V 2|^|K(^) 5o ' a ' 
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From Eq.flMD, via [2 8.146.1] and 8.146.2] and g 0t ± in Eq.fl23]), we get: 



OO OO ; — 

4 - ^ E E E (?) (^M^^-V/) 

\ k=0n=l a,r)=±l x 



2n-l 



T"E E — ^arctanh(p (J (ax / a ) fc 72) if (A.7) 

fe=l <7,??=±1 



„' /IT 'vVago.at i — i i 

Pa ~ V o ' a ~ P+P- 



Therefore this implies: 



oo 1±« 



m a = aJ2 E ^— « fc arctanh(^72) (A.? 

fc=0cr,7?=±l % 

i.e. the second of Eq. (|3.7p . Likewise from the first of Eq. (|2,3p we derive: 

OO 9 _1 

a 



„=±ln=l,*=0 (AQ) 
a=± m=0 ^ «Wr x a 



via [3 8.146.11], hence we obtain the first of Eq. (|3.7p . 

Remark: via identities relating Jacobian elliptic functions it would be possi- 
ble to avoid considering functions that are defined by analytic continuations, 
for instance when the modulus k 2 a > 1 or k'^ < 0, and reduce instead to con- 
sidering only cases with the "simple" arguments (expressing functions by el- 
liptic functions of the arguments h, h' , see for instance the first of Eq. (IA.3j ). 
when necessary): we avoid this as it would hide the nice property that the 
stable and unstable cases are in a sense analytic continuation of each other. 

Summarizing the transformation of coordinates (B,j3)< — >{p' + ,p'_) in 
Eq. (j3.7p is such that t — > p' a e cr ^ 90 ' at generates a solution of the equations 
of motion (if p'± > 0). 

B Polynomials roots 

The proof of the theorem is in the literature: first remark (a classical re- 
sult, see [16l p. 107,301]) that if ao > a% > . . . > a m > 0, and Q m (z) = 
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Y^k=o a k zk ^ 0, then Q m {z) can only vanish if \z\ > 1. Because 

|(1 - z)Q m (z)\ = \a + (ai - a )z + . . . + (a n - a n _i)z n + a n z n+1 \ 

> a - (a - a\)\z\ - ... - (a n _i - a n )2 n - a n \z\ n+1 (B.l) 

> a - (a - ai) - . . . - (a n _i - a n ) - a n = 

Following word by word the proof in p3], write the polynomials P n (z), n > 1 
a + aiz + ...a n 2; n + a n z n+1 + ... + aoz 2n+1 , or (B.2) 
a + aiz + ... a n _iz n_1 + -a n z n + -a n z n + a n _iz n+1 + . . . + a z 2n 

as Q n (z) + 2 n Q* (^r) with n* = n + 1 or n* = n respectively and 

n 

Q n (z) = J2 akz\ Q* n (z)= f z n Q n (z- 1 ) = z n Q n (z- 1 ), (B.3) 

fc=0 

with ao > a\ > . . . > a n (with atf. = for k < n and a n = a n or a n = \a n ). 
Hence Q n ( z ) = aoIK 2 — J") with |^| < 1 by the above remark, and 

Qn(z) = a z n U( l z -4- k ). 

If C is a zero for P n it is Q n (() = -Q n *Q* n {C) or 

icr'nii-c^ 1 ! = icrmc-^i =1 (BA) 
nic-^ 1 ! nii-c**i 

possible only if |£| = 1: because if \a\ < 1 then it is | 1~k z \ > 1 f° r \ z \ > 1 
and < 1 for \z\ < 1 (and n* > 0); or if \z k \ = 1 Eq. (|B.4p means |C| n * = 1. 



C Tree expansion 

To invert the relation x = z C{z) 2 let f(z) d = z C(z) 2 -z = ±±^z 2 + . . . = 
Y^n=2 fn z n , see Eq, (|6.4p . with f n = P n -i a polynomial in r 2 of degree n— 1. 
Consider the equation 

z = x — f(z), with z = x + h(x), or /i(x) = —f(x + /i(x)) (C.l) 

and look for a solution /i(x) = J2T=2 hk% k - 

Let be a tree graph with a root and nodes v into each of which merge 
s v = 0, 1, 2, 3, 4, . . . branches oriented towards the root; each node v carries 
a label n v = 2, . . . with the restriction J2 V s v = k unless v is an end node in 
which case n v = 1 and s v = 0. 
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Such labeled tree graphs will form a family Q(k) of trees identified by 
the rule that two trees are regarded as identical if they can be overlapped by 
pivoting the branches avoiding that any two branches overlap in the course 
of the process. 

Define the "value", Val(0), of 9 as Val(0) = a(6) U vE e ©/n. with 

def 

a (@) = Uvn n >i(~ 1) an d fi = 1- Notice that if he coefficients defining 
f n , n > 2, are all negative then Val(0) > for all 9. Then 

oo 

h(x) = Val (#)^ fc (C- 2 ) 

k=i 6»ee(fc) 

as it is checked by induction (for instance). 

In our case f n — -P n _i(r 2 ) and the contributions of the value of a tree 
9 to the coefficient of x k is a product of symmetric polynomials P ni (r 2 ) = 
Cn') ^h=o Ch) C(n t -h)) r ' 2h w ith symmetrically decreasing coefficients and of 
total degree p — 1 in r 2 ; consequently the polynomials in r 2 being sums of 
products of polynomials enjoying the z< — >z~ 1 symmetry property do enjoy 
it as well. However positivity and monotony is in general lost. 

Notice that Eq. (l6.7j) . (18.2 suggest that the coefficients of x n with n > 
1 are negative polynomials p n {i' 2 )'- in fact —p n (x) are either squares of 
(1 — r 2 ) times polynomials with positive coefficients for n odd or, for n even, 
(1 — r 2 ) 2 times (1 + r 2 )p n (r 2 ) which after performing the multiplication 
become a polynomial with positive coefficients. This property, however, 
does not follow from the above tree expansion and we wonder whether it 
continues to be true at the higher orders. 

Acknowledgment: We are indebted to G. Gentile, A. Giuliani for discus- 
sions and to A. Giuliani for the suggestions quoted in SecEJ 
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